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ABSTRACT 


A  South  China  Sea  warm-core  eddy  with  sea  surface  temperatures  (SST)  higher  than  29.5  "C, 
recently  reported  by  Chu  and  Chang  (1995),  appears  in  the  central  South  China  Sea  (west  of  Luzon 
Island)  in  boreal  spring,  and  strengthens  until  the  onset  of  the  summer  monsoon  (mid-May),  then 
weakens  and  disappears  at  the  end  of  May.  Although  its  size  and  intensity  varies,  the  warm-core 
eddy  releases  large  moisture  and  heat  fluxes  into  the  atmosphere  and  in  turn  affects  the  monsoon 
circulation.  The  transient  features  and  interannual  variabilities  of  the  warm-core  eddy  have  not  yet 
been  studied.  In  this  study  we  use  189,059  temperature  profiles  from  the  U.S.  Navy’s  Master 
Oceanographic  Observation  Data  Set  (MOODS)  to  investigate  the  eddy  features.  We  use  Optimal 
Interpolation  to  create  a  ten-day  interval  synoptic  data  set  for  1964-84  on  a  0.5 °x  1°  grid  (finer 
resolution  in  zonal  direction)  from  the  MOODS  SST  data.  Then,  we  formulate  an  Emprical 
Orthogonal  Functions  (EOF)  analysis.  The  first  EOF  mode  accounts  for  35.5%  of  the  variance  and 
is  characterized  by  a  field  which  is  similar  to  the  mean  field.  The  second  EOF  mode  accounts  for 
21.4%  of  the  variance  and  is  characterized  by  a  meso-scale  eddy  (1 16°-1 19°E,  16°-19°N)  to  the 
west  of  Luzon  Island.  The  size  of  the  eddy  is  around  100,000  km^.  The  time  series  of  EOF2  behave 
differently  between  El  Nino  and  non-EI  Nino  years:  a  warm-core  eddy  often  appears  during  the  non- 
El  Nino  years,  and  a  cold-core  eddy  usually  exists  during  the  El  Nino  years.  A  new  scenario  of  the 
air-sea  feedback  is  proposed  for  the  warm-core  eddy  formation  in  non-El  Nino  years  and  the  cold- 
core  eddy  formation  in  El  Nino  years.  A  numerical  simulation  from  a  three  dimensional  primitive 
equation  model  shows  the  formation  of  a  warm-core  eddy  under  climatological  forcing  (more  likely 
for  the  non-El  Nino  condition). 
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L  INTRODUCTION 


The  South  China  Sea  (SCS)  has  a  bottom  topography  (Figure  1)  that  makes  it  a 
imique  semi-enclosed  ocean  basin  that  is  overlaid  by  a  pronounced  monsoon  surface 
wind.  Extensive  continental  shelves  (less  than  100  m  deep)  are  found  on  the  western  and 
southern  parts,  while  steep  slopes  with  almost  no  shelves  are  found  in  the  eastern  part  of 
SCS.  The  deepest  water  is  confined  to  a  bowl-type  trench.  The  maximum  depth  is  around 
5,000  m. 

Based  on  limited  data  sets,  studies  show  that  eddies  in  the  SCS  are  predominantly 
cyclonic  in  winter  and  spring  and  anti-cyclonic  in  summer,  with  sizes  ranging  from  small 
to  meso-scale  (Huang  et  al.,  1994).  Both  cold  and  warm  eddies  exist  in  the  SCS.  Dale 
(1956)  and  Uda  and  Nakao  (1974)  reported  a  cold  eddy  off  the  central  Vietnamese  Coast 
in  summer.  Nitan  (1970)  foimd  a  cold  eddy  located  to  the  northwest  of  Luzon.  Reports 
from  the  South  China  Sea  Institute  of  Oceanology  (SCSIO,  1985)  indicate  that  in  the 
central  SCS,  a  warm-core  eddy  appears  in  summer  and  winter,  but  more  closer  to 
Vietnam  in  summer  at  the  surface. 

During  boreal  spring,  the  surface  atmospheric  circulation  pattern  changes  from  the 
pattern  of  the  northern  winter  monsoon  to  that  of  the  northern  summer  monsoon.  In 
winter,  a  semi-permanent  surface  subtropical  high  pressure  center  prevails  over  the 
northwestern  Pacific  ivith  its  east-west  axis  located  north  of  20  °N  (Figure  2a).  Cold 
surges  tend  to  move  directly  from  southeastern  China  to  the  SCS  with  north-northeasterly 
winds.  The  winter  cold  surges  tend  to  enhance  the  atmospheric  convective  activity  in  the 
Malaysia-South  China  Sea  Region  (Cheang,  1980). 

During  spring,  the  surface  atmospheric  subtropical  high  pressure  center  moves 
closer  to  the  equatorial  region,  with  its  east- west  axis  between  15°N  and  20  °N  over  the 
northern  region  of  the  SCS  and  northern  Indochina  Peninsula  (Figure  2b).  Cold  surges  do 
not  move  directly  out  from  southeastern  China.  They  travel  a  long  distance  across  South 
Korea  and  southern  Japan  and  around  the  subtropical  high  in  the  western  Pacific.  They 
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arrive  in  the  southern  region  of  the  SCS  as  northeasterlies.  Such  a  semi-permanent 
atmospheric  surface  anticyclone  results  in  downwelling  in  the  upper  ocean. 

Recently  Chu  and  Chang  (1995)  reported  a  warm  pool  (SCSWP)  in  the  central 
SCS  (west  of  Luzon  Island)  during  the  late  spring  season,  that  strengthens  until  the  onset 
of  the  summer  monsoon  (mid-May),  then  weakens  and  disappears  at  the  end  of  May. 
Although  its  size  and  intensity  varies,  the  SCSWP  releases  large  moisture  and  heat  fluxes 
into  the  atmosphere  and  in  turn  affects  the  monsoon  circulation. 

The  thermal  feature  of  eddies  (warm-core  or  cold-core)  greatly  impacts  on  the 
monsoon  circulation.  Therefore,  it  is  very  important  to  know  what  kind  of  eddies  occupy 
the  central  SCS  in  winter  and  spring  seasons .  Can  we  identify  these  eddies  from  the  most 
complete  data  set  -  the  U.S.  Navy’s  Master  Observational  Oceanographic  Data  Set 
(MOODS)?  What  are  the  major  features  of  the  eddies? 

The  data  set  that  we  have  used  is  a  subset  of  MOODS  and  includes  189,059 
profiles  binned  for  the  desired  area  coverage.  The  data  set  covers  the  years  1929  through 
1994  and  comes  from  numerous  instrument  types  including  XBT’s,  CTD’s  and  reversing 
thermometers.  Our  study  field  covers  the  area  5°S  -  25 °N  and  105 °E  -  120°N.  The  data 
were  biimed  for  every  ten  days  creating  756  consecutive  temporal  groups  called  “points”. 
The  Optimum  Interpolation  (01)  technology  was  applied  to  create  a  gridded  data  set  and 
the  GMT  Software  (Wessel  and  Smith,  1993)  was  used  to  plot  SST  contours. 

After  the  gridded  data  set  has  been  established,  we  divided  our  study  into  three 
parts.  First,  a  composite  analysis  of  the  gridded  SST  data  was  done  to  obtain  the 
ensemble  mean  SST  field  and  the  mean  seasonal  variation.  The  latter  shows  strong 
evidence  of  the  existence  of  a  warm  pool  in  the  central  SCS  before  the  summer  monsoon 
onset  (April  and  early  May).  Second,  the  transient  features  of  the  SST  field  were 
quantitatively  investigated  using  the  method  of  Emprical  Orthogonal  Function  (EOF) 
analysis.  An  air-sea  feedback  hypothesis  was  proposed  after  the  EOF  analysis.  Third,  in 
order  to  investigate  the  mechanism  of  the  SCS  warm-core  eddy  formation  during  boreal 
winter  and  summer  under  climatological  forcing  and  its  possible  effects  on  the  summer 
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monsoon  onset,  a  primitive  equation  model  with  second-order  turbulence  closure 
(Blumberg  and  Meller,  1987)  was  used. 

This  thesis  is  organized  as  follows:  Chapter  II  depicts  the  data  sources,  including 
the  U.S.  Navy’s  Master  Oceanographic  Observation  Data  Set  (MOODS)  and  its  temporal 
and  spatial  distributions  as  well  as  the  establishment  of  gridded  data  set.  Chapter  III 
describes  seasonal  and  interannual  variabilities.  Chapter  IV  discusses  numerical 
simulations. 
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Figure  1.  Geography  and  Isobaths  in  Meters,  Showing  the  Bottom  Topography  of 
the  South  China  Sea. 


Figure  2.  Mean  Atmospheric  Surface  Streamline  Analysis  for  (a)  Winter,  and  (b) 
Early  Spring  (From  Cheang,  1980). 
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II.  DATA  SOURCES 


A.  MASTER  OCEANOGRAPfflC  OBSERVATION  DATA  SET  (MOODS) 

The  MOODS  is  a  compilation  of  observed  ocean  data  worldwide  consisting  of  (a) 
temperature-only  profiles;  (b)  both  temperature  and  salinity  profiles,  (c)  sound-speed 
profiles,  and  (d)  surface  temperature  (drifting  buoy).  These  measurements  are,  in  general, 
irregular  in  time  and  space.  In  this  study,  we  analyze  temperature  profiles  obtained  during 
1964-84  from  a  variety  of  instruments.  Due  to  the  shear  size  (more  than  six  million 
profiles)  and  enormous  influx  of  data  to  NAVOCEANO  from  various  sources,  quality 
control  is  a  difficult  task.  Our  study  domain  includes  the  area  5°S  to  25°N  and  105°  to 
120  °E;  the  data  set  within  this  region  consisted  of  189,059  temperature  profiles  after 
rejecting  certain  data  during  quality  control.  The  primary  editing  procedures  included 
removal  of  profiles  with  obviously  erroneous  location,  profiles  with  large  spikes,  and 
profiles  displaying  features  that  do  not  match  the  characteristics  of  surrounding  profiles. 
In  shallow  water,  this  procedure  can  be  partially  automated  but  also  involves  subjective 
interpretation  because  of  the  undersampling  of  MOODS  compared  to  the  spatial  and 
temporal  variability  of  the  oceanography. 

The  temporal  and  spatial  distribution  of  the  MOODS  data  are  irregular.  The 
sparsity  of  profiles  near  the  Sunda  Shelf  and  western  coastal  region  of  Borneo  is  shown  in 
Figure  3  and  the  temporal  distribution  from  1964  to  1984  is  shown  in  Figure  4.  Figure  4 
indicates  that  most  observations  were  obtained  during  the  Vietnamese  War  (1965-69). 
There  are  nearly  14,800  profiles  (maximum  number  of  observations  in  a  year)  in  1968 
(Figure  5).  The  minimum  number  of  observations  is  in  1984  (near  3,000  profiles)  (Figure 
6).  These  figures  indicate  that  the  number  of  observations  within  a  given  year  exhibits  a 
seasonal  variation. 
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Figure  3.  Spatial  Distribution  of  MOODS  Data  during  1964-1984 
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Figure  4.  Temporal  Distribution  of  MOODS  Data  during  1964-1984 


Figure  5.  Temporal  Distribution  of  MOODS  Data  for  1968 
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Figure  6.  Temporal  Distribution  of  MOODS  Data  for  1984 


The  main  limitation  of  the  MOODS  data  is  its  irregular  distribution  in  time  and 
space.  Certain  periods  and  areas  are  over  sampled  while  others  lack  enough  observations 
to  gain  any  meaningful  insight.  Vertical  resolution  and  data  quality  are  also  highly 
variable  depending  much  on  instrument  type  and  sampling  expertise.  There  are  two  data 
sparse  area:  the  Sunda  Shelf  area  and  the  western  coastal  region  of  Borneo  (Figure  3). 
There  are  frequent  ten  to  20  day  gaps  of  no  observations  in  the  whole  South  China  Sea. 
Spatial  and  temporal  irregularities  along  with  the  lack  of  data  in  certain  regions  must  be 
carefully  weighted  in  order  to  avoid  statistically-induced  variability. 

B.  ESTABLISHMENT  OF  GRIDDED  DATA 

Our  near  ten-day  interval  data  set  is  distributed  evenly  in  time  and  unevenly  in 
space.  We  use  the  Optimum  Interpolation  scheme  to  obtain  a  gridded  data  set.  We  use  the 
Pxi'  monthly  SST  climatology  of  Levitus  (1982)  as  a  mean  field.  For  the  same  time 
period,  the  SST  perturbations  T\  ,...,7*^  for  SST  observations  were  obtained  by 
subtracting  the  climatological  value  at  the  closest  grid  point  from  the  observation. 
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optimum  interpolation  assigns  a  weight  to  each  observation  that  accounts  for  variation  in 
spatial  sampling.  The  interpolated  temperature  anomaly  at  the  grid  point,  To,  is  a  linear 
combination  of  the  observed  anomalies,  J*,  jT'j  with  weights 

(1) 

Ml 

The  values  of  the  weights,  a,(/  =  \,2,...,N),  are  found  by  minimizing,  in  a  least  square 
sense,  the  difference  between  the  interpolated  value  at  the  grid  node  and  the  true  value 
there,  and  obtained  from  solving  a  set  of  N  algebraic  equations  (Gardin,  1963) 

N 

=  (2) 

where  z  =  1,2,...,A^ .  Here  X  is  the  signal-to-noise  ratio,  /j.^  is  the  autocorrelation 

between  locations  i  and  j ,  and  is  the  autocorrelation  between  the  grid  node  and 
location  i .  Here,  the  autocorrelation  function  (ACF)  is  defined  by 

=  (3) 

s 

where  denotes  the  independent  space  vector  defining  the  location  of  points  in  a 

sampling  space  R,  and  r  is  the  space  lag  and  the  variance,  /j.  is  computed  by  paring 
the  anomalies  into  bins  depending  upon  their  separation  in  space,  r.  The  values  of  /z  are 
obtained  from  calculating  the  correlation  coefficient  for  all  the  anomaly  pairs  in  each  bin. 
The  ACF  will  be  constructed  for  the  combination  of  different  lags  and  in  the  isotropic 
case 

7(m)  =  7(0)  exp[-./4^  (wAr)^  ]  (4) 
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where  T^rrt)  denotes  the  ACF  value  in  the  bin  with  the  horizontal  separation  /wAr  .  Ar  = 

10  km  is  the  increment  for  the  space  separation.  is  the  horizontal  decorrelation 
scales.  Therefore,  the  autocorrelations  and  are  given  by 


where  niy  and  represent  spatial  separations  between  locations  i  and  j,  and  between 
location  i  and  the  grid  node. 

The  measured  variance  of  the  thermal  fields  is  separated  into  signal  and  noise 
whereby 

5^=5,^+^:  (6) 

The  noise  variance  is  brought  on  from  two  sources,  geophysical  and  instrumentation 
errors.  The  geophysical  error  is  unresolved  thermal  variability  with  scales  smaller  than 
the  typical  time  and  space  scales  between  two  temperature  profiles,  and  the 
instrumentation  error  is  negligible  compared  to  the  geophysical  error  (White  and 
Bernstein,  1979).  In  calculating  the  ACF  value  of  the  first  bin  (zero  to  10  km  lag)  the 
temperature  anomalies  were  never  paired  with  themselves.  The  ACF  value  at  the  zero  lag 
7(0)  was  extrapolated  fi'om  the  first  bin  into  the  origin,  and  therefore  does  not  equal  1 . 
The  signal-to-noise  ratio  is  computed  by  (Sprintall  and  Meyers,  1991). 


^(0) 

1-7(0) 


(7) 


The  larger  the  X ,  the  less  geophysical  error  exists.  If  7(0)  =  1 ,  there  is  no  noise,  =  oo ; 
and  if  7(0)  =  0 ,  there  is  no  signal,  then  A  =  0.  If  /I  >  2 ,  the  ratio  of  the  signal  variance. 


s^. ,  to  the  noise  variance,  s^,  is  greater  than  4,  which  was  considered  a  quite  good 
condition  by  White  et  al.  (1982)  and  Sprintall  and  Meyers  (1991). 

For  each  grid  node  location,  (1)  results  in  a  set  of  N  linear  algebraic  equations  to 
be  solved  for  N  unknowns,  a, ,  by  matrix  inversion.  In  (1)  the  autocorrelations 
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and  the  signal-to-noise  ratio  X  have  been  computed  at  the  surface  from  the  MOODS  data 
set.  Therefore,  the  N  weights  a,  can  be  determined  by  solving  (1),  and  thus  any  new 
observations  can  be  interpolated  into  any  grid  node  G. 

We  applied  Optimal  Interpolation  to  build  up  a  10-day  interval  synoptic  data  set 
for  1964-84  on  a  0.5°  x  1°  grids  (finer  resolution  in  zonal  direction)  from  the  MOODS 
temperature  data.  Our  gridded  data  set  consists  of  756  consecutive  SST  fields  from  1-10 
January  1964  to  21-31  December  1984.  An  example  of  SST  contours  in  1980  is  shown 
in  Figure  7. 
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Figure  7.  Consecutive  SST  Contours  after  the  Establishment  of  Gridded  Data  in 
1980:  (a)  April  1-10,  (b)  April  11-20,  (c)  April  21-30,  (d)  May  1-10,  (e)  May  11-20,  (f) 
May  21-31.  The  Solid  (Dash)  Lines  Means  SST  Higher  (Lower)  than  28°C.  The  Dots 
Indicate  the  MOODS  Stations. 
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in.  SEASONAL  AND  INTERANNUAL  VARIABILITIES 


A.  SCS  WARM-CORE  EDDY  CASE  STUDIES 

Chu  and  Chang  (1995)  reported  a  warm-core  eddy  (SST  ^29 °C)  appears  in  the 
central  SCS  before  the  summer  monsoon  onset  from  analyzing  the  MOODS  data  for 
1966.  A  warm-core  eddy  often  appears  in  spring  and  disappears  in  later  May.  The  SST 
field  in  1966  shows  the  following  features: 

1.  Flat  SST  Pattern  in  March 

This  period  is  right  after  the  cessation  of  the  northern  winter  monsoon  but  some 
hint  of  the  winter  monsoon  characteristics  are  still  present  (Figure  8a-c).  Surface 
isotherms  in  March  show  a  pattern  of  northeast-southwest  oriented  isotherms  with  a 
positive  temperature  gradient  towards  the  southeast  near  the  equator.  The  warmest  water 
(>29°C)  appears  off  the  northwest  coast  of  Borneo.  SST  decreases  from  29°C  near  the 
Borneo  coast  to  24°  at  the  southeast  coast  of  China. 

2.  Early  Stage  of  Warm-Core  Eddy 

This  period  is  in  the  center  of  the  transition  season  (Figure  8d).  In  early  April,  a 
warm  pool  with  the  size  of  near  500,000  km^  appears  near  the  west  coast  of  Luzon.  The 
temperature  of  the  warm  water  increases  from  28  °C  (Figure  8d  and  e)  to  30°C  (Figure 
8f). 

3.  Fully  Developed  Stage  of  SCS  Warm-Core  Eddy 

This  period  occurs  in  late  spring  before  the  summer  monsoon  in  the  SCS  begins 
(Figure  8g).  The  warm  water  pool  expands  southwestward.  During  May  1-10,  the  central 
SCS  was  dominated  by  warm  water  (>30°C).  A  warm  pool  with  SST  higher  than  30.5 °C 
is  located  at  1 1 1  °-l  16°E,  7°-12°N.  This  warm  pool  provides  a  surface  heating  source  for 
the  formation  of  the  atmospheric  low  pressure  center,  which  will  contribute  to  the 
summer  (southwest)  monsoon  onset. 
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Figure  8.  South  China  Sea  SST  Evolution  in  1966:  (a)  March  1-10,  (b)  March  11-20, 
(c)  March  21-31,  (d)  April  1-10  (From  Chu  and  Chang,  1995). 
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Figure  8.  South  China  Sea  SST  Evolution  in  1966:  (e)  April  11-20,  (f)  April  21-30, 
(g)  May  1-10,  (h)  May  11-20  (From  Chu  and  Chang,  1995) . 
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4.  Disappearing  Stage  of  SCS  Warm-Core  Eddy 

This  period  is  just  after  the  summer  monsoon  onset  (Figure  8h).  The  SST  for  May 
11-20  is  relaxed  back  to  the  ensemble  mean  pattern  and  the  SCSWP  disappears.  Warm 
water  (>3 1  °C)  appears  off  the  northwest  coast  of  Borneo  and  decreases  from  Borneo  to 
the  southeast  coast  of  China. 

5.  Cold-Core  Eddy  Building-up  by  Coastal  Upwelling 

Under  the  steady  summer  (southwest)  monsoon  forcing,  cold  water  will  be 
upwelled  near  the  southeast  coast  of  China  and  Vietnam.  A  cold-core  eddy  appears  off 
the  southeast  China  coast  (<25°C)  and  off  the  Vietnamese  coast  (<28°C). 

B.  COMPOSITE  ANALYSIS 

The  existence  of  the  central  SCS  warm-core  eddy  is  likely  to  be  important  for  the 
air-sea  heat  and  moisture  exchanges  and  its  effects  on  the  summer  monsoon  onset.  The 
warm-core  eddy  was  reported  to  the  west  of  Luzon  Island  during  the  winter-to-summer 
monsoon  transition  in  1966  (Chu  and  Chang,  1995).  We  now  examine  the  data  to  see  if 
we  can  obtain  a  warm-core  eddy  signal  from  historical  data.  To  answer  this  question,  we 
need  to  perform  a  composite  analysis.  After  the  gridded  data  is  established,  the  SST  data 
are  represented  by  r(x,,y^.,T^,r,),  where  (X;,y^.)are  the  horizontal  grids,  1964, 
1965,...  1984,  is  the  time  sequence  in  years,  and  t,=\,  2, ...  36,  the  time  sequence  in  ten 
day  intervals.  Before  investigating  the  annual  variation  of  SST,  we  define  the  following 
two  temporal  averages: 

Ar^  1984 -1964  =  21  years  (8) 

which  is  the  annual  mean  value,  and 

T{x,,yj)  (9) 

which  is  the  ensemble  mean  value.  The  ensemble  mean  (1964-1984)  SST  field  over  the 
SCS  (Figure  9)  shows  the  pattern  of  northeast-southwest  oriented  isotherms  with  a 
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positive  temperature  gradient  towards  the  southeast  near  the  equator.  The  ensemble  mean 

has  a  rather  weak  horizontal  temperature  gradient.  T  decreases  from  29  °C  near  the 
Borneo  coast  to  26°C  near  the  southeast  China  coast. 

The  aimual  mean  value  relative  to  the  ensemble  mean, 

T(Xi,yj,t,)=  T(x^,yjJ,)-T(Xi,yj)  ,  shows  the  composite  features  of  the  annual  mean 

SST  anomalies  (Figure  10).  The  annual  mean  SST  anomaly  (T)  has  the  following 
features: 

1.  Cooling  Phase 

In  March,  T  is  negative  almost  everywhere  throughout  the  whole  SCS,  and 
decreases  from  the  west  coast  of  Borneo  to  the  southeast  coast  of  China  (Figure  lOa-c). 
This  indicates  that  the  SCS  is  cooling  during  the  winter  and  early  spring  seasons  and  that 
the  rate  of  cooling  is  weakest  near  the  Borneo  coast  and  strongest  near  the  Taiwan  Strait 
(<-3°C). 

2.  Early  Warming  Phase 

In  early  April  (Figure  lOd),  positive  T  were  found  to  the  west  of  Luzon  Island 
and  the  west  coast  of  Borneo.  At  the  Luzon  Strait  T  is  always  negative  until  late  April 
(21-30  April).  This  means  that  the  early  SCS  warming  starts  in  the  central  SCS  by  some 
local  effects  (e.g.,  a  warm  eddy  but  not  by  the  Kuroshio  intrusion!)  and  at  the  south  SCS. 
It  also  implies  that  the  Kuroshio  water  entering  the  SCS  in  April  might  be  cooler  than  the 
SCS  water  (at  least  at  the  surface). 

3.  Central  SCS  Warming  Phase 

In  May  (Figure  lOg-h),  the  central  SCS  warm-core  eddy  becomes  stronger  than 
the  early  warming  phase.  During  11-20  May  (Figure  lOh),  the  central  SCS  warm-core 
eddy  (  T  >  1 .8  °C)  accounts  for  the  maximum  annual  mean  warm  anomaly. 
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Figure  9.  The  Ensemble  Mean  SST  Field  (1964-84).  The  Solid  (Dash)  Lines  Means 
SST  Higher  (Lower)  than  28°  C. 
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Figure  10.  Average  SST  Anomalies  Relative  to  the  Ensemble  Mean  for  (a)  March  1- 
10,  (b)  March  11-20,  (c)  March  21-31,  (d)  April  1-10,  (e)  April  11-20,  (f)  April  21-30. 
Solid(Dash)  Contours  Denote  Positive(Negative)  Anomalies. 
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Figure  10.  Average  SST  Anomalies  Relative  to  the  Ensemble  Mean  for  (g)  May  1-10, 
(h)  May  11-20,  (i)  May  21-31,  G)  June  1-10,  (k)  June  11-20,  (1)  June  21-30. 
Solid(Dash)  Contours  Denote  Positive(Negative)  Anomalies. 
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4.  Mature  Warming  Phase 

After  the  summer  monsoon  onset,  the  June  T  becomes  positive  almost  everywhere 
throughout  the  whole  SCS,  and  increases  from  the  west  coast  of  Borneo  (a!0.6°C)  to  the 
southeast  coast  of  China  («2.4°C).  This  indicates  that  after  the  summer  monsoon  onset, 
the  SCS  warming  is  weakest  near  the  Borneo  coast  and  strongest  near  the  Taiwan  Strait 
(>2.4°C). 


C.  EMPRICAL  ORTHOGONAL  FUNCTION  (EOF)  ANALYSIS 


The  composite  analysis  demonstrated  the  existence  of  a  central  SCS  warm-core 
eddy  in  the  annual  mean  SST  anomaly  data.  What  are  the  transient  features  of  this  meso- 
scale  eddy?  What  are  the  synoptic  and  interannual  variabilities  of  the  eddy?  We  can  use 
the  Emprical  Orthogonal  Function  (EOF)  technique  to  answer  these  questions. 

1.  EOFs  for  the  SST  Anomalies 
The  SST  synoptic  anomalies  obtained  by 


are  re-arranged  into  a  NxP  matrix,  7’(r„,^),n=  1,2,...,Y;  and  />=  1,2,...,P .  Here  P  = 

756  is  the  total  number  of  time  points  used  for  computing  the  covariance  matrix,  i.e.,  21 
years  of  every  10-day  measurements;  N  =  255,  corresponds  to  the  number  of  grids 
Q  =  1,2,.  ..,17;/  =  1,2,.  ..15) .  Lorenz  (1956)  used  the  eigenvectors  of  the  covariance  matrix 
associated  with  a  sample  data  field  to  calculate  the  time  series  of  coefficients  or 
predictors  for  the  statistical  weather  prediction.  This  method,  called  EOF  analysis, 
separates  the  data  sets  into  eigenmodes.  Generally  speaking,  each  mode  has  an  associated 
variance,  a  non-dimensional  spatial  pattern,  and  a  dimensional  time  series.  From  this  data 
matrix  a  255-square  spatial  covariance  matrix  is  calculated  by 
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where  n  and  m(\,2,...,25S)  denote  the  grid  locations.  The  diagonal  elements  of  the 
covariance  matrix  ={n  =  \,2,...,N)  are  the  variance  at  location  r„.  The  off-diagonal 
elements  are  the  covariance  with  spatial  lag  equal  to  the  difference  between  the  row  and 
column  indices.  This  symmetric  matrix  has  its  255  real  eigenvalues  ,  and  eigenvectors 
(paQ'j)  such  that 

=  (  =  1A...A'  (12) 

The  eigenvectors  called  Emprical  Orthogonal  Functions.  Each  <f>^  is  a 

255-point  (17x15  grid)  distribution  of  SST  anomaly  pattern.  The  eigenvalues, 
X^{a  =  \,2,...,N),  are  all  positive  and  the  summation  of  them,  ,  equals  the  total 
variance.  Therefore,  is  considered  as  the  portion  of  total  variance  “explained”  by  the 
EOF  .  It  is  convenient  to  label  the  eigenfunctions  (j>^  so  that  the  eigenvalues  are  in 
descending  order,  i.e., 

A,  >2.2  >^3  >...  (13) 

2.  Principal  EOF  Modes 

Intraseasonal  and  interannual  anomalies  can  be  quantitatively  investigated  with 
the  method  of  emprical  orthogonal  function  (EOF)  analysis  and  a  number  of  its 
generalized  forms.  In  order  to  delineate  the  major  modes  of  variability  in  the  winter  and 
spring  seasons  (December  to  May)  of  21  years  (1964-84),  we  perform  a  conventional 
EOF  analysis  and  obtain  the  first  six  leading  EOFs.  The  first  six  leading  EOF  modes  are 
able  to  accoimt  for  81.7%  of  the  total  variance  during  this  season.  They  are  summarized 
in  Table  1.  When  combined,  EOFl  and  EOF2  explain  over  half  (57%)  of  the  total 
variance.  Each  EOF  mode  is  normalized  so  that  its  total  spatial  variance  is  equal  to  unity. 
So,  those  patterns  of  the  first  six  EOFs  are  enough  to  explain  the  spatial  anomalies  of 
whole  SCS.  Hence,  it  suffices  to  focus  on  the  first  six  EOFs.  Their  structures  are  shovm 
in  Figure  1 1 .  Notice  that  the  normalized  EOF  values  are  multiplied  by  10. 
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EOF 

Variance(%) 

Accumulative  Variance 

(%) 

1 

35.6 

35.6 

2 

21.4 

57.0 

3 

12.4 

69.4 

4 

5.3 

74.7 

5 

4.7 

79.4 

6 

2.3 

1 

81.7 

Table  1.  Variances  of  the  First  Six  Leading  EOFs. 


The  EOFl  mode  has  an  orderly  homogeneous  pattern  of  correlation  and  accounts 
for  up  to  35.6%  of  the  total  spatial  variance  during  the  winter  and  spring  seasons 
(December  to  May).  The  NE  -  SW  oriented  isotherm  pattern  over  this  region  is  believed 
to  be  related  to  the  ensemble  mean  pattern.  The  EOF2  (Figure  11b)  mode  shows  a  strong 
eddy  pattern  with  a  closed  isoline  of  -4°C  (actual  value  is  -0.4  °C)  to  the  west  of  Luzon 
Island  (116°-119°E,  16°-19°N)  and  accounts  for  about  21.4%  of  variance.  The  size  of 
the  eddy  is  around  100,000  km^.  However,  we  still  cannot  identify  whether  the  eddy  in 
EOF2  represents  a  warm  or  cold  ring  since  the  values  in  Figure  8  are  normalized.  We 
must  therefore  examine  the  products  of  the  EOF  values  and  the  corresponding  time  series 
of  the  amplitudes. 
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Figure  11.  The  First  Six  EOF  Modes  for  the  SST  Anomalies  (°C)  Computed  From 
December-May  During  1964-84;  (a)  EOFl,  (b)  EOF2,  (c)  EOF3,  (d)  EOF4,  (e) 
EOF5,  (f)  EOF6.  The  EOF  Model  Values  have  been  Multiplied  by  10. 


3.  Temporal  Variabilities 

The  data  matrix,  T{r„,Tp)  is  approximately  written  by 

=  (14) 

where  B^(tp)  is  a  set  of  time  series  of  size  P,  representing  the  temporal  variation  of  the 
associated  spatial  pattern  described  by  EOF  ■  The  result  of  the  time  series  seems  to 
be  a  projection  of  the  SST  anomaly  through  a  ‘filter’  of  EOF  modes  during  a  time  scale. 

The  time  series  of  EOFl  amplitudes  B^(7p)  for  1964-84  are  shown  in  Figure  12. 
The  EOFl  mode  is  always  negative  throughout  the  whole  South  China  Sea 

(Figure  11a).  Therefore,  5,(^)>0,  corresponds  to  negative  SST  anomalies,  and 
5,(?^)<0,  corresponds  to  positive  SST  anomalies.  Taking  the  year  of  1968  as  an 
example,  5,(?^)  is  always  positive  from  early  January  to  early  May,  which  implies 

negative  SST  anomalies  in  the  whole  SCS  during  this  period.  Furthermore,  the  strength 
of  the  negative  SST  anomalies  is  weak  near  the  west  coast  of  Borneo  and  increases 
northwestward  towards  the  southeast  coast  of  China  (Figure  1  la). 

The  time  series  of  EOF2  amplitudes  B2(Tp)  for  1964-84  are  shown  in  Figure  13. 
The  EOF2  mode  (Figure  11b)  is  almost  negative  throughout  the  whole  SCS 

except  near  Hainan  Island  and  shows  a  strong  eddy  with  closed  isoline  of  -4°C  to  the 
west  of  Luzon  Island.  As  before,  ^2  (^)  >  0,  corresponds  to  negative  SST  anomalies  (or  a 
cold-core  eddy),  and  B2(tp)<0,  corresponds  to  positive  SST  anomalies  (or  a  warm-core 
eddy).  Taking  the  year  of  1964  as  an  example,  B2{7p)  is  almost  always  negative  from 
early  January  to  early  May  and  reaches  a  minimum  value  (-15)  in  1-10  April,  which 
implies  positive  SST  anomalies  over  the  whole  SCS  and  a  warm-core  eddy  to  the  west  of 
Luzon  Island  during  the  1964  winter  and  spring  seasons.  The  strength  of  the  warm-core 
eddy  increases  from  February  to  April.  The  SST  anomaly  of  the  warm-core  eddy  in  April 
is  around  6°C  (-0.4°C  x  -15).  The  actual  SST  anomaly  (6°C)  is  the  product  of  the 
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second  EOF  mode  (-0.4  °C)  and  its  corresponding  amplitude  (-15).  On  the  other  hand, 
taking  the  year  of  1972  as  another  example,  BjiJp)  is  almost  always  positive  from  early 

January  to  early  May  with  a  mean  value  around  3,  which  implies  negative  SST  anomalies 
over  the  whole  SCS  and  a  cold-core  eddy  to  the  west  of  Luzon  Island  during  1972  winter 
and  spring  seasons.  The  strength  of  the  cold-core  eddy  fluctuates  from  January  to  May. 
The  SST  anomaly  of  the  cold-core  eddy  in  1972  is  around  -1.2°C  (-0.4°C  x  3). 

D.  SCS  COLD-CORE  EDDY  DURING  EL  Nmo  YEARS 

From  the  above  analysis,  we  see  a  different  thermal  behavior  in  the  central  SCS 
between  the  1964  and  1972  winter  monsoon  seasons:  a  warm-core  eddy  in  1964  and  a 
cold-core  eddy  in  1972.  What  might  cause  this  interannual  variability?  Initially,  one 
might  think  of  a  change  in  the  large-scale  surface  atmospheric  circulation.  A  reasonable 
guess  is  the  El  Nino  event  (1972  is  an  El  Nino  year).  During  an  El  Nino  year,  the 
monsoon  circulation  was  greatly  weakened  (Kuo,  1989).  The  wind  stress  curl  over  the 
SCS  might  be  different  between  non-El  Nino  and  El  Nino  years.  Following  this  line  of 
reasoning,  we  averaged  the  amplitudes  of  EOF  1  and  EOF2  for  El  Nino  years  and  non-El 
Nino  years,  respectively.  The  set  of  El  Nino  years  consists  of  1965,  1969,  1972,  1976, 
1982,  and  1983,  and  the  set  of  non-El  Nino  years  includes  the  rest  of  the  years.  The  mean 
fields  for  non-El  Nino  years  generally  may  represent  conditions  closer  to  the 
climatological  mean.  The  averaged  amplitudes  of  EOFl  of  non-El  Nino  years  indicate 
positive  values  from  February  to  March  (Figure  14a),  which  implies  negative  SST 
anomalies  over  all  the  SCS  during  this  period.  Furthermore,  the  strength  of  the  negative 
SST  anomalies  is  weak  near  the  west  coast  of  Borneo  and  increases  northwestward 
towards  the  southeast  coast  of  China.  Conversely,  the  averaged  amplitudes  of  EOFl  of  El 
Nino  years  are  negative  values  (Figure  14b).  Its  SST  anomaly  pattern  is  opposite  to  the 
pattern  of  the  non-El  Nino  years  during  this  period  of  time. 

The  averaged  amplitudes  of  EOF2  of  non-El  Nino  years  indicate  negative  values 
during  the  whole  winter  and  spring  seasons  (Figure  15a),  which  implies  positive  SST 
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anomalies  over  all  the  SCS  and  a  warm-core  eddy  to  the  west  of  Luzon  Island  during  this 
period.  However,  the  averaged  amplitudes  of  EOF2  during  El  Nino  years  indicate 
positive  values  during  the  whole  winter  and  spring  seasons  (Figure  15b),  which  implies 
negative  SST  anomalies  over  all  the  SCS  and  a  cold-core  eddy  to  the  west  of  Luzon 
Island  during  this  period.  The  averaged  amplitude  is  usually  smaller  than  individual 
years.  This  implies  that  a  SCS  warm-core  eddy  appears  frequently  from  January  to  May 
during  non-El  Nino  years. 

Om  data  set  also  shows  that  during  El  Nino  years  a  cold-core  eddy  will  frequently 
appear  from  January  to  May  in  the  central  SCS  (Figure  15b).  After  relating  the  difference 
in  the  SST  anomaly  pattern  to  non-El  Nino  and  El  Nino  years,  we  believe  that  the 
appearance  of  a  warm-core  or  a  cold-core  eddy  is  related  to  the  global  atmospheric 
circulation. 

E.  SCENARIOS  FOR  THE  CENTRAL  SCS  EDDY  FORMATION 

The  MOODS  data  set  shows  that  during  boreal  winter  and  spring  a  warm-core 
eddy  often  appears  in  the  central  SCS  except  for  El  Nino  years  when  a  cold-core  eddy 
usually  appears.  To  examine  this  further,  we  use  climatological  data  sets  such  as  the 
monthly  mean  surface  wind  stress  data  (Hellerman  and  Rosenstein,  1983)  and  monthly 
mean  temperature  and  salinity  data  sets  (Levitus,  1982)  to  study  the  eddy  formation 
mechanisms  in  non-El  Nino  years  since  the  mean  features  of  the  non-El  Nino  years  are 
more  likely  closer  to  the  climatology. 
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Figure  12.  Time  Series  of  EOFl  Amplitudes  (January  to  May)  for  1976-84. 
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Figure  13.  Time  Series  of  EOF2  Amplitudes  (January  to  May)  for  1964-75. 
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1977  EOF2 
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Figure  13.  Time  Series  of  E0F2  Amplitudes  (Januaiy  to  May)  for  1976-84. 


Non-El  Nino  years  average  EOF1 


El  Nino  years  average  EOF1 


Figure  14.  Time  Series  (January  to  May)  of  EOFl  Amplitudes  Averaged  Anomaly: 
(a)  Non-Nino  Years,  (b)  and  El  Nino  Years  (1965, 1969, 1972, 1976, 1982, 1983). 


34 


Non-Elnio  years  average  EOF2 


Elnio  years  average  EOF2 


Figure  15.  Time  Series  (January  to  May)  of  EOF2  Amplitudes  Averaged  Anomaly^ 
(a)  Non-Nino  Years,  (b)  and  El  Nino  Years  (1965, 1969, 1972, 1976, 1982, 1983). 
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We  try  to  answer  the  following  questions:  What  is  the  mechanism  for  a  warm- 
core  eddy  formation  under  climatological  conditions  (more  likely  representing  non-El 
Nino  years)?  What  is  the  difference  of  the  surface  forcing  over  the  SCS  between  non-El 
Nino  and  El  Nino  years?  What  is  the  mechanism  for  a  cold-core  eddy  formation  in  El 
Nino  years? 

Two  possible  mechanisms  may  explain  the  central  SCS  warm-core  eddy 
formation  during  the  winter  monsoon  season:  (1)  the  effects  of  downwelling,  and  (2)  the 
Kuroshio  intrusion.  Figure  lOd  indicates  that  the  early  SCS  warming  starts  at  the  central 
and  south  SCS,  not  at  the  Luzon  Strait.  This  implies  that  the  Kurashio  intrusion  is  not  a 
key  factor  to  cause  the  central  SCS  warm-core  eddy  formation.  Therefore,  the 
downwelling  effect  is  the  predominant  mechanism  for  the  central  SCSWP  formation  in 
the  boreal  winter  and  spring. 

1.  Ekman  Pumping  in  Central  SCS 

The  monthly  climatological  wind  data  (Hellerman  and  Rosenstein,  1983)  over  the 
SCS  and  adjacent  areas  were  used  to  compute  the  Ekman  pumping  velocity  (Figures  16- 
19). 


f 


(15) 


Here  are  the  surface  wind  stresses.  Figures  16-19  shows  the  Ekman  pumping 

velocity  by  the  color  contours  for  March-Jime.  The  arrows  indicate  the  surface  wind 
stress.  We  see  that  during  March-May  (Figure  16-18)  a  downwelling  region  wdth  a 
vertical  velocity  around  -0.025  cm/s  appears  in  the  central  SCS.  After  the  summer 
monsoon  onset  (June),  the  surface  winds  change  direction,  and  a  weak  upwelling  area 
with  vertical  velocity  around  0.005  cm/s  appears  in  the  central  SCS  between  14°N  and 
1 9  °N  (Figure  19). 
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2.  Thermal  Feature  of  the  Central  SCS 

Figiire  20  shows  two  sets  of  April  MOODS  temperature  profiles.  One  is  for  1968 
(a  non-El  Nino  year),  and  the  other  is  for  1982  (an  El  Nino  year).  The  central  SCS  water 
is  always  highly  stratified  before  the  summer  monsoon  onset. 

The  ocean  mixed  layer  temperature  is  much  higher  than  the  sublayer  temperature. 
The  Ekman  downwelling  in  the  central  SCS  during  the  winter  and  spring  seasons  will 
suppress  the  vertical  cold  advection  from  the  sublayer.  This  downwelling  effect,  plus  the 
increase  in  solar  radiation  in  spring  (April-May),  may  provide  some  important  pre¬ 
conditions  for  the  central  SCS  warm-core  eddy  formation. 

3.  An  Air-Sea  Feedback  for  a  Warm-Core  Eddy  (Non-El  Nino) 

We  propose  an  air-sea  feedback  scenario  for  the  central  SCS  warm-core  eddy 
formation  in  non-El  Nino  years  (Figure  21).  From  late  winter  to  spring,  an  anticyclone 
appears  in  the  central  SCS  from  the  ensemble  mean  atmospheric  surface  streamline 
analysis  (Figure  22a).  This  anticyclone  generates  downwelling  in  the  central  SCS  and  in 
turn  prevents  the  cold  deep  water  from  being  advected  to  the  surface.  This  will  promote 
the  formation  of  the  warm-core  eddy  in  the  central  SCS.  As  the  central  SCS  warm-core 
eddy  persists,  being  around  1  °C  warmer  than  the  surroundings,  an  atmospheric  surface 
low  pressure  center  will  be  generated  above  the  warm  water,  which  may  trigger  the 
summer  monsoon  onset. 

After  the  summer  monsoon  onset,  an  atmospheric  surface  cyclone  will  occupy  the 
central  SCS  (Figure  22b).  This  cyclone  will  generate  upwelling  in  the  central  SCS,  which 
entrains  the  deep  cold  water  into  the  surface  mixed  layer.  This  upwelling  effect  will 
finally  destroy  the  central  SCS  warm-core  eddy. 

During  El  Nino  events,  the  Asian  monsoon  circulation  is  greatly  weakened  (Kuo, 
1989)  and  so  will  be  the  downwelling  effect.  Therefore,  the  warm-core  eddy  is  not  likely 
to  occur  in  the  boreal  winter  and  spring  during  El  Nino  events.  The  mechanism  for  a 
cold-core  eddy  formation  in  El  Nino  years  should  be  further  studied. 
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Figure  16.  The  Ekman  Pumping  Velocity  During  March.  The  Colors  Associated 
With  Positive  (Negative)  Values  Indicate  Upwelling  (Downwelling)  (cm/s).  The 
Arrows  Denote  the  Mean  Wind  Stresses  (From  Hellerman  and  Rosenstein,  1983). 
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Figure  17.  The  Ekman  Pumping  Velocity  During  April.  The  Colors  Associated  With 
Positive  (Negative)  Values  Indicate  Upwelling  (Downwelling)  (cm/s).  The  Arrows 
Denote  the  Mean  Wind  Stresses  (From  Heilerman  and  Rosenstein,  1983). 
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Figure  18.  The  Ekman  Pumping  Velocity  During  May.  The  Colors  Associated  With 
Positive  (Negative)  Values  Indicate  Upwelling  (Downwelling)  (cm/s).  The  Arrows 
Denote  the  Mean  Wind  Stresses  (From  Hellerman  and  Rosenstein,  1983). 
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Figure  19.  The  Ekman  Pumping  Velocity  During  June.  The  Colors  Associated  With 
Positive  (Negative)  Values  Indicate  Upwelling  (Downwelling)  (cm/s).  The  Arrows 
Denote  the  Mean  Wind  Stresses  (From  Hellerman  and  Rosenstein,  1983). 
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Figure  22.  Mean  Atmospheric  Surface  Streamline  Analysis  for  (a)  Early  Spring,  and 
(b)  Summer  (From  Cheang,  1980). 
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IV.  NUMERICAL  SIMULATIONS 


We  use  a  numerical  model  to  simulate  the  formation  of  the  central  SCS  warm- 
core  eddy.  Through  this  simulation,  we  hope  to  verify  the  downwelling  effect  scenario. 

A.  PRIMITIVE  EQUATION  MODEL 

The  model  we  use  is  the  three  dimensional  model  developed  by  Blumberg  and 
Mellor  (1983,  1987)  with  the  hydrostatic  and  Boussinesq  approximations  (Bryan,  1969) 
and  has  the  following  features:  (1)  horizontal  curvilinear  coordinates  and  an  "Arakawa  C" 
scheme  (Arakawa  and  Lamb,  1977),  (2)  sigma  coordinates  in  the  vertical  with  realistic 
bathymetry,  (3)  a  free  surface,  (4)  a  second-order  turbulence  closure  model  for  the 
vertical  viscosity  (Mellor  and  Yamada,  1974, 1982),  (5)  horizontal  diffusivity  coefficients 
calculated  by  the  Samagorinsky  (1963)  parameterization,  and  (6)  split  time  steps  for 
barotropic  and  baroclinic  modes.  The  basic  equations  are  given  by 

du  du  ^  \  ^  d  du. 

=  (16) 

dv  ^  s  \  dp  d  dv.  _ 

— -t-UoVv  +  w - fu  = - —  +  — (iCw — )  +  F„ 

a  dz  ^  ^  dz^  ^  dz^  ^ 


pg  =  — 


VoF  +  — =  0 
& 

dd  de  d de^  ^ 

—  +  F  oV^-i-w —  =  — — )  +  Ff, 
a  oz  oz  oz 

a  a  a  ”  a  ^ 

where  V  =  (w,v)  is  the  horizontal  velocity,  w  the  vertical  velocity,  and  g  the  gravitational 
acceleration,  is  the  characteristic  density,  p  the  in  situ  density,  0  the  temperature,  S 
the  salinity, /?  the  pressure,  the  vertical  eddy  diffusivities  for  turbulent  mixing 
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of  momentum  {K^^)  and  temperature  and  salinity  {K„),  (F^,Fy,Fg,Fs)  the  horizontal 
mixing  terms  representing  the  subgrid  scale  processes.  The  density  p  is  computed  by  the 
equation  of  state  of  the  form 

P  =  P(0,S)  (17) 

given  by  Fofonoff  (1962). 

The  model  domain  is  99°E-121°E,  3°S-25°N  which  includes  the  South  China 
Sea  and  the  Gulf  of  Thailand.  With  its  20  km  horizontal  resolution  and  23  vertical  sigma 
coordinate  levels,  it  has  125x162x23  grid  points  (Figure  23).  All  depths  less  than  10  m  in 
the  domain  are  set  to  be  10  m  in  order  to  satisfy  the  constraint  77+  77  >  0 .  In  the  vertical 

the  domain  is  divided  into  23  levels  <T„(n  =  l,2,...,23).  Here  cr„  =nAa,A(T  =  ~j^^.  In 

the  simulation,  we  use  the  Samagorinsky  parameterization  to  compute  the  horizontal 
mixing  coefficients  'with  C=0.1,  which  results  in  the  horizontal  viscosity  ranging  from 
200  to  500  mV^  The  background  vertical  mixing  coefficient  is  lO^m^s**  (Blumberg  and 
Mellor,  1983,  1987). 

B.  LATERAL  BOUNDARY  CONDITIONS 

There  are  several  major  straits  connecting  the  SCS  to  the  Pacific  Ocean  and 
surrounding  seas.  Many  are  very  shallow  (i.e.,  the  Strait  of  Malacca)  or  irregular  (i.e.,  the 
Balabac  Strait)  and  difficult  to  be  adequately  represented  in  a  numerical  model. 
Therefore,  we  close  the  Strait  of  Malacca,  and  all  the  small  straits  between  Luzon  Island 
and  Borneo  in  the  numerical  model.  We  combine  the  Karimata  Strait  and  the  Gasper 
Strait  into  one  open  boundary.  After  this  treatment,  our  model  has  three  open  boimdaries: 
a  northern  boimdary  (Taiwan  Strait),  an  eastern  boimdary  (Luzon  Strait),  and  a  southern 
boundary  (Karimata-Gasper  Strait).  The  transports  at  the  three  open  boimdaries  were 
determined  according  to  Wyrtki's  (1961)  estimation  with  some  modification  due  to  the 
boundary  treatment.  Table  2  shows  the  seasonal  variation  of  the  transport 
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(li'v-lO  ms  )  at  the  three  open  boundaries.  The  summation  of  the  total  transport 
through  the  open  boundaries  is  always  zero  in  order  to  satisfy  mass  conservation. 


LONGITUDE 


Figure  23.  Model  Domain  of  South  China  Sea  with  20  km  x  20  km  Grids. 
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Month 

Feb. 

Aug. 

Oct. 

Dec. 

Gasper  &  Karimata  Straits 
(Eastward  Positive) 

0.0 

-4.0 

-3.0 

1.0 

4.3 

Luzon  Strait 
(Eastward  Positive) 

-3.5 

0.0 

2.5 

-0.6 

■ 

Formosa  Strait 
(Northwards  Positive) 

-0.9 

0.0 

1.0 

0.5 

-0.4 

-0.9 

Table  2.  The  Bi-monthly  Variation  of  Mass  Transport  (Sv)  at  The  Open  Boundaries. 
The  Values  were  Taken  from  Wyrtld  (1961). 
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The  barotropic  velocity  at  the  open  boundaries  is  obtained  by 

(18) 

where  /  is  the  boundary  width,  and  n  means  the  normal  direction  to  the  boundary.  The 
baroclinic  velocity,  Vf^(a)  is  calculated  by 

=  (19) 

1-e 

During  the  numerical  integration,  both  and  were  interpolated  into  the 
time  step  from  the  bi-monthly  values  of  the  volume  transport  listed  in  Table  2. 

On  the  lateral  open  boundaries  temperature  (T)  and  salinity  (S)  are  prescribed  at 
the  inflow  boundaries,  whereas  at  the  outflow  boundaries,  the  advection  equation, 
namely, 

^(7-,S)+F,^(7’,S)  =  0  (20) 

a  on 

is  solved  for  both  barotropical  and  baroclinic  modes. 

C.  STABILITY  CONSTRAINTS 

The  stability  constraints  were  derived  from  various  linearized  subsets  of  equations 
(Blumberg  and  Mellor,  1987).  Since  zlx=29  km,  Hjfiax=5000  m,  and  Ufnax  ~  1  the 
typical  CFL  condition  derived  from  Blumberg  and  Mellor  (1987)  for  the  external  mode  is 
45  s  and  for  the  internal  mode  is  2000  s.  In  our  numerical  model,  the  time  step  was  set  as 
25  s  for  the  external  mode  and  as  900  s  for  the  internal  mode. 

D.  SURFACE  FORCING 

The  atmospheric  forcing  includes  wind  forcing  and  thermodynamic  forcing.  We 
took  the  monthly  mean  climatological  wind  stress  data  (Hellerman  and  Rosenstein,  1983) 
as  the  wind  data  for  the  middle  of  the  month  and  interpolated  them  for  each  day.  The 
wind  stress  has  a  typical  magnitude  of  1  dyne/cm  (corresponding  to  a  surface  wind  speed 


of  10  m/s).  Northeasterly  winds  dominate  during  the  winter  and  spring  seasons  while 
southwesterly  winds  prevail  during  the  summer  monsoon  season. 

There  are  two  approaches  to  setting  up  the  surface  thermal  forcing;  (1)  flux 
forcing  when  the  surface  heat  and  salt  fluxes,  m‘^ )  and  ^^(m  )  are  given,  and 

(2)  restoring  forcing  when  the  fluxes  are  not  given  and  the  model  surface  temperature  (T) 
and  salinity  (S)  are  relaxed  to  the  observed  values  (Fob)  and  (Sob)  '• 

^TP 

“  ^\Qh  (21) 

oz 

K„^  =  a,Q,*afi(S„-S). 

OZ 

where  Kuvs.  the  thermal  exchange  coefficient  and  C  is  the  relaxation  constant  taken  to  be 
0.7  m/d,  which  is  equivalent  to  a  relaxation  time  of  43  days  for  an  upper  layer  30  m  thick. 
{a^,a^  are  (0,l)-type  switcher  parameters:  when  a,  =  haj  =  0,  only  flux  forcing  is 
present;  a,  =  0,^2  =  X  only  the  restoring  forcing  is  present,  or,  =  Totj  =  h  both  flux 
forcing  and  the  restoring  forcing  are  present.  Process  studies  usually  use  only  the 
restoring  forcing  as  the  surface  thermal  forcing  (e.g.,  Brady  and  Gent,  1994).  Therefore, 
in  this  simulation  we  use  only  the  restoring  forcing.  We  construct  the  data 

through  interpolation  firom  the  monthly  mean  temperature  and  salinity  data  (Levitus, 

1982).  Both  the  horizontal  and  vertical  resolutions  of  the  Levitus  (T,S)  data  are  different 
from  our  model  resolution.  Hence,  we  use  optimal  interpolation  (01)  to  convert  the 
Levitus  (T,S)  data  to  our  model  grids. 

E.  MODEL  SIMULATED  SST  FIELDS 

The  model  was  initiated  from  April  mean  (T,S)  fields  (Levitus,  1982)  and  forced 
by  the  surface  wind  stress  (Hellerman  and  Rosenstein,  1983)  and  thermodynamical  fluxes 
depicted  in  (21).  The  model  was  integrated  for  three  years.  The  third  year’s  outputs  were 
used  for  discussion. 

Three  days  of  the  simulated  surface  currents  were  picked  for  investigating  the 
time  evolution  of  the  circulation  pattern.  We  choose  30  January  for  the  winter  monsoon 
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forcing,  30  March  for  the  boreal  spring  forcing,  and  30  April  for  the  pre-summer 
monsoon  forcing. 

On  30  January  a  strong  southward  coastal  jet  is  present  near  the  Vietnam  coast, 
and  no  cyclonic  or  anticyclonic  rotation  is  observed  in  the  central  SCS  (Figure  24a).  The 
isolines  of  SST  in  the  central  SCS  are  almost  parallel  to  the  southeast  China  coast  (Figure 
25a). 

On  30  March  the  strength  of  the  southward  Vietnam  coastal  current  is  greatly 
reduced.  An  anticyclone  with  a  strong  northward  branch  is  generated  between  13°N- 
18°N,  111  °E-1 15°E  (Figure  24b).  Several  small  cyclones  are  also  simulated,  e.g.,  a  small 
and  weak  one  (next  to  the  anticyclone)  is  located  at  13  °E-16°E,  1 15°E-1 16°E.  A  warm- 
core  eddy  is  also  simulated  associate  with  the  anticyclone  (Figure  25b). 

On  30  April  the  southward  coastal  jet  disappears  near  the  Vietnam  coast.  In  the 
SCS  an  anticyclone  persists  between  13°N-18°N,  111°E-115°E  (Figure  24c).  A  warm- 
core  eddy  is  also  simulated  associate  with  the  anticyclone  (Figure  25c). 

The  simulations  confirm  our  scenario  for  the  warm-core  eddy  formation  under  the 
climatological  forcing.  The  cold-core  eddy  formation  during  El  Nino  events  should  be 
studied  in  the  future  research. 
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Figure  24.  Model  Simulation  of  Surface  Circulation  for  (a)  30  January,  (b)  30 
March,  and  (c)  30  April. 


Figure  25.  Model  Simulation  of  SST  (°C)  for  (a)  30  January,  (b)  30  March,  and  (c) 
30  April. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  goal  of  our  study  was  to  detect  the  meso-scale  eddies  in  the  central  South 
China  Sea  during  the  winter  and  spring  seasons  from  the  Navy’s  Master  Observational 
Oceanographic  Data  Set  (MOODS),  and  to  investigate  mechanisms  for  the  eddy 
formation.  Optimal  interpolation,  a  composite  analysis  and  an  EOF  analysis  were  used  to 
study  the  first  objective;  a  3-D  numerical  model  was  used  for  the  second  objective. 

After  quality  control  procedures  were  employed,  the  data  set  was  composed  of 
180,059  observed  profiles  spanning  the  years  1964  to  1984.  Due  to  a  sparsity  of  data,  ten- 
day  bins  were  constructed  for  each  year.  Optimal  inteipolation  (01)  was  used  to  establish 
a  gridded  SST  data  set  for  the  SCS.  We  applied  optimal  interpolation  to  construct  a  ten- 
day  interval  synoptic  data  set  (756  total)  from  1-10  January  1964  to  21-31  December 
1984  on  a  0.5°  X  1°  grid  (finer  resolution  in  the  zonal  direction)  from  the  MOODS 
temperature  data.  The  gridded  data  size  is  45x25x756. 

The  composite  analysis  on  the  gridded  data  set  indicates  that  the  annual  mean  SST 
seasonal  variation  has  four  stages:  (1)  a  cooling  phase,  (2)  an  early  warming  phase,  (3)  a 
central  SCS  warming  phase,  and  (4)  a  mature  warming  phase.  A  warm-core  eddy  starts  to 
occur  in  the  early  warming  phase  and  persists  imtil  the  end  of  the  central  SCS  warming 
phase.  It  disappears  in  the  mature  warming  phase. 

The  EOF  analysis  of  the  SST  anomaly  data  shows  an  eddy  in  the  central  SCS 
being  present  in  the  second  principal  mode.  The  eddy  has  distinctive  thermal  features 
between  non-El  Nino  and  El  Nino  years.  During  non-El  Nino  years  the  central  SCS  has  a 
warm-core  eddy.  However,  during  El  Nino  years  the  central  SCS  has  a  cold-core  eddy. 
This  implies  that  the  eddy  formation  in  the  central  SCS  is  related  to  the  global 
atmospheric  circulation. 
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A  hypothesis  for  the  SCS  eddy  formation  under  the  climatological  forcing  (more 
likely  non-El  Nino  years)  was  also  proposed  in  this  study.  During  late  winter  an 
atmospheric  anticyclone  over  the  central  SCS  generates  downwelling  and  in  turn  prevents 
the  cold  deep  water  from  being  advected  to  the  surface.  As  solar  radiation  increased  in  the 
late  winter  monsoon,. the  central  SCS  surface  water  is  heated  more  than  the  surrounding 
waters  by  the  downwelling  effect.  This  will  promote  the  formation  of  a  warm-core  eddy 
in  the  central  SCS.  As  the  central  SCS  warm-core  eddy  persists,  being  a  few  degrees 
warmer  than  its  surroundings,  an  atmospheric  surface  low  pressure  center  will  be 
generated  above  the  warm  water,  which  may  trigger  the  summer  monsoon  onset.  After 
the  summer  monsoon  onset,  an  atmospheric  surface  cyclone  occupies  the  central  SCS. 
This  cyclone  will  generate  upwelling  in  the  central  SCS,  which  advects  the  deep  cold 
water  into  the  surface  mixed  layer.  This  upwelling  effect  will  finally  destroy  the  central 
SCS  warm-core  eddy.  This  dovrawelling  hypothesis  on  the  SCS  warm-core  eddy 
formation  during  non-El  Nino  years  was  verified  by  a  3-D  primitive  equation  model.  The 
air-sea  feedback  scenario  has  not  been  verified  for  El  Nino  events. 

B.  RECOMMENDATIONS 

We  recommend  that  further  studies  should  incorporate  (1)  3-D  thermal  features  of 
the  SCS  meso-scale  eddies,  (2)  constructing  a  composite  for  the  SCS  surface  wind  fields, 
to  test  if  they  are  upwelling  or  downwelling  favorable  wind  field;  (3)  running  the  3-D 
numerical  model  xmder  the  El  Nino  vdnd  forcing. 
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APPENDIX  A.  TIME  SERIES  OF  EOF  3  AMPLITUDES  FOR  1967-84 


1964  EOF3  1965  EOF3  1966  EOF3 
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1970  EQF3  1971  EOF3  1972  EQF3 


1973  EOF3  1974EOF3  1975  EOF3 


1976  E0F3 


1977  E0F3 


1978  E0F3 


1979  E0F3  1980  E0F3  1981  EOF3 


1982  EOF3  1983  EOF3  1984  EQF3 
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APPENDIX  B.  TIME  SERIES  OF  EOF  4  AMPLITUDES  FOR  1967-84 
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APPENDIX  D.  TIME  SERIES  OF  EOF  3  AMPLITUDES  AVERAGE  ANOMALY 


El  Nino  years  average  EOF3 


Non-Ei  Nino  years  average  EOF3 
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APPENDIX  E.  TIME  SERIES  OF  EOF  4  AMPLITUDES  AVERAGE  ANOMALY 


El  Nino  years  average  EOF4 


Non-El  Nino  years  average  EOF4 
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APPENDIX  F.  TIME  SERIES  OF  EOF  5  AMPLITUDES  AVERAGE  ANOMALY 

El  Nino  years  average  EOF5 


Months(10  days/pt) 


Non-El  Nino  years  average  EOF5 


Months(10  days/pt) 
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